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ABSTRACT 

We report a detection of water vapor in the protoplanetary disk around DoAr 44 with the Texas 
Echelon Cross Echelle Spectrograph — a visitor instrument on the Gemini north telescope. The DoAr 
44 disk consists of an optically thick inner ring and outer disk, separated by a dust-cleared 36 AU gap, 
and has therefore been termed “pre-transitional”. To date, this is the only disk with a large inner gap 
known to harbor detectable quantities of warm (T ^ 450 K) water vapor. In this work, we detect and 
spectrally resolve three mid-infrared pure rotational emission lines of water vapor from this source, 
and use the shapes of the emission lines to constrain the location of the water vapor. We find that the 
emission originates near 0.3 AU — the inner disk region. This characteristic region coincides with that 
inferred for both optically thick and thin thermal infrared dust emission, as well as rovibrational CO 
emission. The presence of water in the dust-depleted region implies substantial columns of hydrogen 
(> lO^^cm”^) as the water vapor would otherwise be destroyed by photodissociation. Combined 
with the dust modeling, this column implies a gas/small-dust ratio in the optically thin dusty region 
of > 1000. These results demonstrate that DoAr 44 has maintained similar physical and chemical 
conditions to classical protoplanetary disks in its terrestrial-planet forming regions, in spite of having 
formed a large gap. 

Subject headings: stars: pre-main sequence — protoplanetary disks 


1. INTRODUCTION 

Detecting protoplanets in their parent disks would 
allow us to directly study the interactions and rela¬ 
tionships between planets and their birth environments. 
A sub-class of protoplanetary disks termed transition 
disks remain strong candidates for disks that harbor gi¬ 
ant (proto)planets, and thus have generated a signifi¬ 
cant amount of study. Transition disks were originally 
identified via a dearth of near-IR excess emission, in- 
dica tive of depletion of small grains close to the star 


fe-Fi- 

Strom et alJll989l:IKoerner et alJll993HCalvet et alJ 

|20Q2 

), and a leading explanation for the formation of 


many of these dust-depleted disks is clearing by tidal 
interactions with giant planets fe.g. iRice et al.l 120031 
iZhu et ani2011l : iDodson-Robinson fc Salvkll2011ll . Many 
transition disks have since been imaged at millime¬ 


ter wavelengths, confi rming the presence of inner dust 
holes or clearings le.g. lBrown et al.ll2009l [Andrews et al.l 
1201 If) . In addition, there have been several detections 


of companions in transition disks using 
masking and/or direct imaging (e.g. 

sparse aperture 

Kraus & Ireland 

2012 

:lHuelamo et alJl2011l:IBiller et al.ll2012l:IClose et al. 

2014 

)■ 


Transition disks are far from homogeneous, and their 
inner regions can have complex properties. Near- 
IR photometry and interferometry show that many 
transition disks have dust in their inner regions, in¬ 
cluding some, terme d pre-transitionah with optically 
thick inner rims (e.g. lEisner et al.ll200^ lEsnaillat et al.l 
[Moll . Many transition disks also have azimuthally 
asymmetric dust and gas structures (iMuto et al.ll20I2t 

Ivan der Marel et ^120131 iCasassus et al.ll20I3[ l. suggest- 
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ing that they are actively being sculpted by planets. But 
assessing the potential for planet growth, and character¬ 
izing planets’ birth environments, also requires knowl¬ 
edge of the gas content and chemistry. 

A powerful probe of conditions in the inner regions 
of transition disks is water vapor. Water vapor is eas¬ 
ily photodissociated by FUV radiation from young, ac¬ 
creting stars. Therefore, the presence of water va¬ 
por requires shielding of UV radiation and/or suffi¬ 
ciently high gas densities for chemical production (e.g. 
lAdamkovics et all 12014 iDn fc BerginI 1201^. with pos¬ 
sible modifications due to mixing (e.s. iCiesla fc Cu^ 
l2006t[Mbertsson et al.ll201^ . Spitzer-IRS spectra of nine 
transition disks (reduced according to the procedure de¬ 
scribed in iPontoppidan et al.ll2010H) . a water-rich classi¬ 
cal disk, and a water vapor model, are shown in Figure 
[H demonstrating that nearly all transition disks show 
no evidence for warm water vapor. Given that ^40% 
of young low-mass stars with disk s show water vapor 
emission (jPontoppidan et aT1l2010bll . it is unlikely that 
these transition disks are collectively water-poor through 
chance alone. Instead, the clearing of the inner regions of 
transition disks likely reduces dust-induced UV shielding 
and gas densities which , combined, “dry out” the inner 
disk (| Zhang et al.ll2013ll . In fact, the Spitzer-IRS spec¬ 
trum from TW Hya does show water vapor emission at 
longer wavelengths, but this is emitted by cooler water 
vapor located at the inner edge of the optically thick 
outer d isk, where gas and dust densities retu rn to higher 
values (|Zhang et ai.ll201.‘ll : lKamp et al.ll201.1ll . Therefore, 
water must already be incorporated into planetesimals, 
or be delivered, if water-rich planets are to be produced 
in the inner disk. 

An interesting exception to this trend is the disk 
around the star DoAr 44 (also called ROX 44 and Haro 
1-16). This d i sk wa s classified as pre-transitional by 
lEspaillat et al.l (120101) . with a spectral energy distribu¬ 
tion and near-IR spectrum consistent with the presence 
of an optically thick ring near 0.25 AU, optically thin 
dust out to 2 AU, and a cleared gap between 2 and 36 
AU. Its Spitzer-IRS spectrum shows clear evidence for 
water vapor (Figure [ H. with a H 9 O/CO line ratio sim¬ 
ilar to classical disks (iSalvk et all 12011^ . suggesting a 
similar C/0 ratio in the emitting region. Given the high 
upper level energies (up to 6000 K) for the observed wa¬ 
ter lines, they likely arise from gas close to the young star 
— a hypothesis that can be confirmed using spectrally 
resolved observations. The presence of water vapor in 
this disk presents us with an opportunity to better un¬ 
derstand the conditions in the inner regions of a disk that 
has likely already formed giant planets, and what impli¬ 
cations this may have for still-forming terrestrial planets. 

In this work, we present spectrally resolved observa¬ 
tions of water vapor from DoAr 44 and discuss their im¬ 
plications. In Section [21 we discuss the observations and 
the detection of water vapor. In Section [31 we estimate 
the radial origin of the emission, the water vapor tem¬ 
perature and column density, and the implied hydrogen 
column density. In Section IH we present our conclusions 
and discuss their implications. 

2. OBSERVATIONS AND DATA REDUCTION 

Our target is the protoplanetary disk around DoAr 
44 — a young K3 star in the p Ophiuchus cloud 


([Bouvier fc Appenzelled [199^ . Spectra were obtained 
with the Texas Echelon Cross Echelle Spectrograph 
(TEXES; iLacv et ^ 120021) — a visitor instrument on 
Gemini north — on 2014 August 15 and 16 UT (Pro¬ 
gram ID: GN-2014B-Q-40). The total target integration 
times were 56 and 39 minutes, respectively, on the two 
nights. Observations were performed in the “hi-medium” 
cross-dispersed mode, with a resolution of R^80,000 at 
a wavelength of ~12.4 /rm. We measure a FWHM of 
3.5±0.1kms“^ (R~85000) for a narrow sky absorption 
feature, confirming that we achieve the expected reso¬ 
lution. The slit size was 0.52" x 4", enabling nodding 
along the slit. 

Spectra were reduced utilizing the TEXES pipeline. 
Reduction includes flat fielding, subtraction of nod pairs 
and optimal extraction of spectral traces. The pipeline 
derives the wavelength solution from a single atmospheric 
line in the setting combined with knowledge of the in¬ 
strument optics, and solutions are found to be accurate 
to ~ 0.5kms“^. To correct for telluric absorption, the 
bright asteroid Eunomia was observed at a similar air- 
mass, and the spectra from DoAr 44 were divided by 
this spectrum. For absolute flux calibration, we scale 
to the 0.55 Jy 12.4 / im coii tinuum flux level measured by 
iPontoppidan et aP (|2010bD in each order. Uncertainty in 
absolute flux calibra tion for Spitzer-IRS spectra is ^10% 
([Furlan et al.l l200^ : in addition, transition disks often 
show v ariability at the 10%- 20% level in this wavelength 
region (jEspaillat et al.ll201lD . Therefore, there is a 10- 
20% uncertainty in absolute flux. 

The full wavelength coverage was 12.36-12.46 p,m. Fig¬ 
ure m shows the portion of the fully reduced spectrum 
containing three water vapor emission lines expected 
to emit strongly. All three water vapor lines are de¬ 
tected, with line/continuum ratios near 0.2. Details are 
shown in Table 111 Because the emission lines have low 
line/continuum ratios and are broad, we have binned 
our spectrum to 2x10“'’’ pm (5kms“^) — a factor of 
5 larger than the native pixel size of ^4xl0“^/rm 
(0.9kms“^). Data are missing in gaps between orders. 

Vertical lines mark expected line locations. Emission 
line centers are well separated from telluric water, but the 
wide line profiles from DoAr 44 result in some overlap. 
The atmospheric emission here causes a lower signal-to- 
noise ratio (S/N) on the blue side of the line profiles, 
which we believe accounts for the observed line asymme¬ 
try. Alternatively, the disk emission lines could be in¬ 
trinsically asymmetric, with strongly redshifted emission 
or blueshifted absorption. However, this source is accret¬ 
ing a t a moderate rate (9 x 10“® Mr7^vr~^: lEspaillat et al.l 
I 2 OIOD and has no known evidence for an outflow. In ad¬ 
dition, the broad emission lines are similar in shape to 
observed GO rovibr ational emission lines from this source 
(jSalvk et al.l 1201 fbh . which show no evidence for asym¬ 
metry. Therefore, we assume here that the lines arise 
from the protoplanetary disk itself and that any asym¬ 
metries are insignificant. 

3. RESULTS 

3.1. Location of water vapor emission 

In FigurejSl we show a weighted (by S/N) average wa¬ 
ter vapor line profile, as well as an averaged CO rovibra¬ 
tional emission line profile (|Salvk et al.ll20lfbl) . While 
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Wavelength [fim] 


Fig. 1.— Portions of continuum-subtracted 5pit2er-IRS spectra from transition disks, a water-rich disk around DR Tau (divided by 3 
for clarity) and an LTE slab model (T=500 K, A^H 20 = 4 x 10^® cm~^ and R = 0.6 AU). Transition disk classifications from fCalvet et ^ 
II2002I) ; IFurlan et al.l ll2006li ; IBrown et al.l ll2007li ; lEspaillat et al.l ||2010I~) . 



Wavelength [//m] 


Fig. 2.— Observed spectrum (gray), binned to a resolution of 2 A, and atmospheric transmission (green). Dashed vertical lines mark 
expected location of water vapor emission. Labels show upper and lower state quantum numbers (J/CaKc) state energies. A red 

dotted line marks the 0.56 Jy continuum. 
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TABLE 1 


Line 

(J/C„kU 

(/im) 

Eupper 

(K) 

Spin 

Flux'" 

I 64 13 —>■ 15l 14 

12.3765 

4948 

ortho 

9.25 ± 0.94 

174 13 —>■ I 63 14 

12.3971 

5781 

ortho 

3.99 ± 0.94 

I 63 13 —>■ 152 14 

12.4079 

4945 

para 

3.26 ± 0.96 


^ Theore tical line cente rs from HITRAN molecular 
database (IRothman et al.1 1201311 . shifted by the Earth- 
induced Doppler shift, ~ 28kms“^, plus the heliocentric 
velocity of the source, — 8 kms~^, measured from rovibra- 
tional CO emission lines llSalvk et al.ll2011bll . 

10 ~^^ erg cm~^ s~^ from Gaussian fit. Error prescription 
from iLenz & Ayres! ll 19921 1. Errors do not include a 10%- 
20 % uncertainty in absolute flux. 

there are differences between the two line profiles, the 
relatively low S/N of the water emission precludes as¬ 
cribing this to true differences in the emitting location of 
the molecules. Instead, we focus on the primary proper¬ 
ties of the emission line — namely, that it is broad, and 
possibly double-peaked. The width of the emission line is 
clearly inconsistent with arising from the outer disk be¬ 
yond 36 AU, as the maximum observed velocity in such a 
case would be only 2.5 km s“^. Thr oughout this work, we 
assume a stellar mass of 1.4 Mq (lAndrews et al.l[2009fl 
and a disk inclination of 25°, based on an elliptical f it to 
the dust emission map in Ivan Dishoeck et al.l ()2015h . 

The peak emission and half width at zero iten- 
sity (HWZI) have corresponding velocities of 20 and 
65kms“^. Assuming a 10kms“^ uncertainty, these im¬ 
ply emission arising from an annulus with inner ra¬ 
dius 0.05to'o2 A-U and outer r adius 0.6tn;^ AU. Follow¬ 
ing the procedure outlined in iSalvk etahl (l2011bD . we 
also compare the observed lineshape to simplified disk 
models in which the line luminosity, L, is modeled as 

L = /^°“* L{R)dR. We assume L{R) oc ® from i?rin 

to i?mid and L{R) oc R~^ from Rmid to Rout, such that 
the emission drops rapidly beyond We take Rout 

as 100 AU. 

Figure [3] demonstrates that a model with i?i„ = 0.05 
AU and i?mid = 0.6 AU is consistent with the data. It 
also shows how the model varies with Rin and Rmid, and 
conclusively demonstrates that the observed lineshape is 
inconsistent with very large values of Rin- Therefore, the 
water vapor emission arises in the inner disk. 

Models in which the emission falls more rapidly with 
radius (i.e., with a steeper L(R) profile) produce the bulk 
of their flux near i?i„, and the observed double-peaked 
profile could be produced by such a model with slightly 
larger Rin- A model of this type, with an inner ra dius of 
0.45 AU, is utilized bv iPontopoMan et al.l ()2ni0all to fit 
water vapor emission from the circumstellar disk around 
the classical T Tauri star RNO 90, and also provides a 
reasonably good fit to the observed emission from DoAr 
44 (see Figure [3]). Correcting for RNO 90’s higher incli¬ 
nation of 37°, the implied inner radius for DoAr 44 would 
be - 0.22 AU. 

In summary, we find that the bulk of the line luminos¬ 
ity emerges near 0.3 AU, either from an annulus extend¬ 
ing from ~0.05-0.6 AU, or from a region with an inner 
radius of 0.22 AU. In the former case, the inner radius 
is close to the ^0.04 AU sublimation radius (assuming 


a sublimation tempe rature of 1500 K and L = l.dL©; 
lAndrews et al.l[2009f) . The two alternative models may 
be distinguished with higher S/N data by the depth of 
the line profile’s central depression. 

3.2. Water temperature and column, and implications 
for local conditions 

We investigate the properties of the water vapor col¬ 
umn by considering it as a single-temperature slab of gas 
in LTE, with a column density, A^H 2 0 : area ttR? (where 
we refer to R as the emitting radius), and temperature, T 
(and assuming a distance of 125 pc and turbulent velocity 
of 2kms“^). The TEXES emission lines alone are not 
constraining, so we co nsider them in con junction with 
the Spitzer-\RS data (ISalvk et al.l I2011an over the 14- 
17 pLvn range. The combined fit requires T > 500 and 
A^HaO ^ 10 ^® cm“^, with the lowest temperatures requir¬ 
ing higher A^HaO and vice-versa. A representative model 
has A^HaO = 10^9 cm-2, T = 800K and R = 0.15 AU 
(the “hot” model in Figure S]). 

A fit to the full Spitzer-WXS spectral range found 
T = 450K, = 2 x 10i®cm-2, and R = 0.79 AU 

(ISalvk et al.ll2011^ . This “warm” model (Figure 14]) un¬ 
derpredicts the TEXES line fluxes by a factor of 10. Con¬ 
versely, the single-temperature model fitting the TEXES 
data does not produce enough long wavelength flux to 
match the Spitzer data beyond ~17 pLui. This discrep¬ 
ancy suggests that a range of temperatures is required 
to explain the observed water vapor emission. 

The modeled column density, A^UoO 10^® cm~^ and 
implied H 2 O/CO ratio of >1 (|Saivk et al.l 12011^ are 
consis tent with those found for classical prot oplanetary 
disks (iCarr fc Naiitall 20 lll : ISalvk et'aI1l2011^ . This col¬ 
umn of water implies shielding from photodissociating 
FUV radiation. Whether this shielding is provided pri¬ 
marily by dust grains or by self-shielding depends on 
the dust opacity, with self-shielding dominating when 
FUV opacities drop bel ow 1% of canonical ISM values 
(iBethell fc Bergiiil l2009fl . According to the models of 
lEspaillat et al.l (|2010f) . the water vapor emitting region 
overlaps both an optically thick dust ring (with radius 
< 0.4 AU) and an optically thin dusty region within 2 
AU, whose vertical optical depth is only 0.036 at 10 /rm. 
For 0.1 /rm pure silicate dus t grain s with optical con¬ 
stants from lOssenkopf et al.l (Il992| l. the corresponding 
optical depth at a photodissociating wavelength of 110 
nm is 0.56. The actual optical depth will depend on 
the precise composition and dust grain size. However, if 
grains have grown beyond 0.1 /rm in size, the UV opacity 
will be lower. Therefore, the dust is unlikely to provide 
sufficient opacity to shield the H 2 O in this region. 

If the dust optical dep th is low, models of self-shielding 
(|Bethell &: BerginI l2009f) relate observed water column 
density to vertical hydrogen column density. For an 
example disk around AA Tau, with an accret ion rate 
of ~ 3 X lO“®M 0 yr“^ (|Gullbrinef et al.l flOOSH , a wa¬ 
ter column density of 10 ^® cm“^ (the minimum col¬ 
umn density required to fit the combined TEXES and 
Spitzer spectra) corresponds to a vertical hydrogen col¬ 
umn of ^ 10^^ cm“^, for a range of temperatures (T > 
600 K). DoAr 44’s slightly higher accretion rate of 9 x 
lO“®M 0 yr“^ (and correspondingly higher accretion lu¬ 
minosity) implies a factor of a few larger hydrogen col- 
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umn. In addition, if th e water vapor emissio n lines are 
sub-thermally excited (|Meiierink et alJl2009f) . required 
water vapor column densities and hence hydrogen col¬ 
umn densities will be higher still. Therefore, the water 
vapor detection implies that it arises from a region with 
a vertical hydrogen column TVh > 10^^cm“^. These 
models assume a canonical water abundance, but true 
abundances may differ from this value by factors of 5- 
10, depend ing on relative rates o f gas and icy body radial 
transport (|Ciesla fc Cuzzil [200611 . If planets form DoAr 
44’s gap, these may block the transport of icy bodies, 
resulting in a lower water abundance and higher implied 
hydrogen column. 

If the water vapor emitting region overlaps with the 
optically thin region p roducing the 10 ^m dust emission 
(jEspailiat et al.lf2010ll . we can estimate the gas/(small)- 
dust ratio in this region. For 0.1-1 /im sized sili- 
cate grains, the opac ity at 10 ^m is ~ 2100cm^g“^ 
(|Ossenkopf et al.lll992ll . so an optical depth of 0.036 im¬ 
plies a dust mass surface density of 1.7 x 10“^gcm“^. 
Dividing the hydrogen mass surface density by the dust 
mass surface density gives a gas/dust ratio of ~1000. In 
this estimate, all of the (optically thin) dust is accounted 
for, but the hydro gen column derived from the models of 
IBethell &: Bergiiil is that lying above the water-emitting 
layer, and so represents a lower limit to the total vertical 
hydrogen column. Therefore, the gas/(small-)dust ratio 
in the optically thin dust region is >1000. 

4. CONCLUSIONS AND DISCUSSION 

We detect velocity-resolved water vapor emission lines 
from the pre-transition disk DoAr 44, using the TEXES 
spectrograph on Gemini north. We show that: 

• The water vapor emission arises from the inner disk 
near 0.3 AU — well inside the ~36 AU inner rim of 
the optically thick outer disk, and overlapping the 
dusty inner disk regions modeled bv lEspaillat et all 

dMa). 

• The water vapor observed with TEXES requires 
warmer temperatures than are required to fit the 
Spitzer-IRS spectrum. Therefore, the emission 
arises from a region with a range of temperatures. 

• The best-fit water column density (A^HaO > 
10^®cm“^) implies a vertical hydrogen column of 
IVh > 10^ ^cm~^. based on mode ls for water self¬ 
shielding (|Bethell fc BerginI 12009(1. For the opti¬ 
cally thin dusty region (|Espaillat et al.ll2010ll . this 
implies a gas/(small-)dust ratio > 1000. 

Many transitio n disks are active accretors (e.g. 
INaiita et al.ll2007ll and harbor CO gas in their inner dust- 

depleted regions (IRettig et ‘^12004 iSalvk etaI1l2007ll — 

facts that support a planet-cleari ng origin for these disks 
(Dodson- Robinson fc Salvkll2011h . Here, we furthermore 
demonstrate that the hydrogen column in the inner disk 
of DoAr 44 must be > 10^^ cm“^ in order to maintain 
a significant column of water vapor. We also derive a 
super-canonical gas/(small-)dust ratio, which may be the 
result of grain growth — a process likely requ ired to ex¬ 
plain the low opacity throughout the disk gap (iZhu et al.l 


Conversely, this work implies that the majority of tran¬ 
sition disks, which do not have detectable warm water 
vapor emission, have not been able to maintain high in¬ 
ner disk gas columns. Note that we would not expect to 
detect water vapor as strong as is observed from DoAr 44 
without also detecting emission with the Spitzer-IRS (for 
spectra with similar S/N to those in Figured]). There¬ 
fore, if other transition disks harbor water vapor in their 
inner regions, the amount of water vapor must be signif¬ 
icantly less than that observed from DoAr 44. 

Why has DoAr 44 maintained a gas-rich, water-rich 
inner region, while other transition disks have not? Al¬ 
though DoAr 44 has been termed pre-transitional, im¬ 
plying that it may be dissipating its inner disk to evolve 
into a transition disk, the expected timescale for dissipa¬ 
tion of the inner disk is short. The viscous timescale at 
0.3 AU for an actively accr eting disk like DoAr 44 (with 
a = 0.01) is ^1000 years (IRudenlfiQQQd — less than a 
hundredth of the disk lifetime. More likely, the inner 
disk is being replenished by the outer disk, in agreement 
with recent observations showing that t he outer regions 
of transition disk gaps are not gas-free (ICasassus et al.l 
120131 Ivan der Marel et al.ll201^ . If DoAr 44 is not evolv¬ 
ing into a transition disk, the differences between DoAr 
44 and other transition disks may lie in the architec¬ 
tures of the planetary systems they harbor. The sim¬ 
ilarity in stellar and disk properties between DoAr 44 
and other transition disks (lAndrews et al.l[20Tlll suggests 
that the late-stage chemistry of the terrestrial planet 
forming region may depend on a stochastic aspect of 
planet growth/disk evolution. 

Continued study of DoAr 44 will shed more light on 
the formation scenario for transition disks, and on im¬ 
plications for inner disk chemistry. Higher S/N mid-IR 
spectroscopy, ideally combined with spectro-astrometry, 
would place firmer constraints on the radial origin of the 
water vapor emission, allowing a uniquely powerful view 
of the terrestrial planet forming region. Complementary 
observations of CO gas in the lO’s of AU region have 
already been obtained with the Atacama Large Millime¬ 
ter Array (van der Marel et al. 2015, in preparation). 
By combining these sets of observations, we will better 
understand the overall architecture of the dust-depleted 
gap, and what this implies for any embedded planets. 
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Fig. 3. — First p anel: Weighted average water vapor line profile (black), binned to 5kms“^, and rovibrational CO line profile (gray; 
ISalvk et al.ll2011bl) . Vertical lines mark line peak and line wing (HWZI) velocities, corresponding to outer and inner radii, respectively. 
Second panel: Weighted average water vapor line profile (black) and emission models with different values of The models with 

= 0.05 and 0.5 AU have Rmid = 0-6 AU; the model with Rin = 36 AU has Rmid = 40 AU. Third panel: Same as second panel, but 
models have R^ = 0.05 and d ifferent values of Rm\d- Fo urth panel: Weighted average water vapor line profile (black) and water vapor 
emission line model (red) from [Pontoppidan et all (I2010al ) with Rin = 0.45 AU. 



Fig. 4. — Top left: Continuum-sub tracted TEXES spec trum (black), “hot” model fit to TEXES and 14-17 pm Spitzerd&ta. (blue), “warm” 
model fit to Spitzer-IKS data from ISalvk et ahl l)2011al i (red) and sum of models (gray). Top right: Continuum-subtracted Spitzer-IKS 
spectrum (black) and combined hot and warm models (gray). Bottom: Continuum-subtracted Spitzer-IRS spectrum (black), hot model 
(blue) and warm model (red). 
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